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ABSTRACT 

The galactic field's late-type stellar single and binary population is calculated on the 
observationally well-constrained supposition that all stars form as binaries with in- 
variant properties in discrete star formation events. A recently developed tool (Marks, 
Kroupa & Oh) is used to evolve the binary star distributions in star clusters for a 
few Myr until an equilibrium situation is achieved which has a particular mixture of 
single and binary stars. On cluster dissolution the population enters the galactic field 
with these characteristics. The different contributions of single stars and binaries from 
individual star clusters which are selected from a power-law embedded star cluster 
mass function are then added up. This gives rise to integrated galactic field binary dis- 
tribution functions (IGBDFs) resembling a galactic field's stellar content (Dynamical 
Population Synthesis). It is found that the binary proportion in the galactic field of 
a galaxy is larger the lower the minimum cluster mass, Meci,min, the lower the star 
formation rate, SFR, the steeper the embedded star cluster mass function (described 
by index P) and the larger the typical size of forming star clusters in the considered 
galaxy. In particular, period-, mass-ratio- and eccentricity IGBDFs for the Milky Way 
are modelled using Meci,min — 5 Mq, SFR= 3 Mq yr~^ and (3 — 2 which are jus- 
tified by observations. For rh « 0.1 — 0.3 pc, the half-mass radius of an embedded 
cluster, the afore mentioned theoretical IGBDFs agree with independently observed 
distributions, suggesting that the individual discrete star formation events in the MW 
generally formed compact star clusters. Of all late-type binaries, 50 per cent stem from 
-^eci ^ 300AfQ clusters, while 50 per cent of all single stars were born in Mod > lO^M© 
clusters. Comparison of the G-dwarf and M-dwarf binary population indicates that the 
stars formed in mass-segregated clusters. In particular it is pointed out that although 
in the present model all M-dwarfs are born in binary systems, in the Milky Way's 
Galactic field the majority ends up being single stars. This work predicts that today's 
binary frequency in elliptical galaxies is lower than in spiral and in dwarf-galaxies. 
The period and mass-ratio distributions in these galaxies are explicitly predicted. 

Key words: Galaxy: stellar content - binaries: general - star clusters: general - open 
clusters and associations: general - solar neighbourhood - methods:numcrical 



1 INTRODUCTION 



The Milky-Way's (MW) Galactic field stellar late- 
type population (spectral-types G to M) consists of 
roughly 50 per cent single stars and 50 per cent 
binaries JPuguennov fc Mavoil Il99ll: Mayor etal] 1 19921 : 



iFischer fc Marcvlll992l : iRaghavan et al.N2010l '). Throughout 



this work we will use the lDuguennov fc Mavoil (|l99ll . here- 
after DM91) results for G-dwarfs as the canonical Galactic 
field population to which data is being compared, since it 
is the only existing long-term (13 yrs) survey performed by 
a single team, i.e. the DM91 data provides a homogeneous 
data set. 
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In a spe ctroscopic survey of halo, thick and thin disc 
populations, ICarnev et al.l (|2005l ) find for high proper mo- 
tion stars that 28±3 per cent of all metal-poor ([Fc/H]^ — 1) 
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centre-of-mass (cm-)system^ are binaries. A similar number 
(26 ± 3 per cent) is found for the more metal-rich stars in 
their sample ([Fe/H]> —1). These values compare well with 
a binary-fraction of 22 per cent identified by DM91 over the 
same period range (1.9 < P < 7500 d). So there is no evi- 
dence for variation of binary properties in the Galactic field 
over cosmological time. 

The origin of the Galactic field composition is a result 
of star formation but has not been predicted with success. 
Indeed, iFished { 20041 ) finds in his theory of isolated star 
formation a bell-shaped binary period distribution function 
(period BDF) broadly similar to that of the Galactic field, 
but is not able to make specific predictions concer ning the 
form o f the period- and o t her B DFs. As discussed in lKroupal 
l|201lD . iMoeckel fc Bat3 l|2010t )'s seminal SPH simulations 
result in a semi-major axis BDF with too many binaries 
with orbits around a few astromical units (AU). Further- 
more, their mass-ratio distribution shows too few binaries 
with q = m-xlm\ < 0.8, where mi and m2 are the primary- 
and secondary-component mass, respectively. The reason for 
this might be a collapse which is too deep resulting in a too 
dense cluster and thus too efficient binary disruption. Ad- 
ditionaly after 10 Myr of dynamical evolution including in- 
stantaneous residual-gas expulsion, orbits with semi-major 
axes g > 1 AU are under-represented. Stellar feedback (e.g. 
lBatdl2009l ). i.e. self-regulated star formation, might help to 
remedy this by reducing the depth of the collapse through 
an increased pressure. 

Indeed, direct cloud collapse calculations are very lim- 
ited in predicting binary star properties owing to the se- 
vere computational difficulties of treating the magneto- 
hydrodynamics together with correct radiation transfer and 
evolving atomic and molecular opacities during collapse. 
Such computations are also prohibitively expensive. There 
is therefore currently no sufficient numerical framework to 
derive the multiplicity properties of stars. While such com- 
putations are important for understanding the physical pro- 
cesses in observed star formation regions, they do not have 
predicitive power yet. In particular, larger systems, e.g. 
galaxies, cannot be synthesized. 

Furthermore, the outcome of star formation computa- 
tions cannot be expected to result in the binary properties of 
the Galactic field since clustered star fo rmation is the dom- 
inan t discrete star formation event (e.g. lLada fc Ladall2003l : 
iLada 2010; Bressert et al. 201Q) and binaries are dynami- 
cally processed in star clusters before they become part of 
a galactic field. Therefore, a galactic field stellar population 
is the sum over all discrete star formation events, which on 
dissolution contribute a number of single stars and binaries 
depende nd on the star formation cond itions and dynamical 
history (|Kroupalll995bl : iGoodwinllioiol ) . This is the topic of 
the present investigation. 

Observations of low-density pre-main sequence pop- 
ulations, i.e. "distributed" star formation, show that 
most, if not all stars form as members of binaries 
and that they exhi bit a period BDF that is rising to- 
wards long periods JSimon et al.l Il995l: iGhez et aL _ 19971: 



20081 : iKroupa fc Petr-Gotzensll201lD . Indeed, there is a sim- 



Kohler fc Leinerd Il998l: iDuchend Il999l : IConnellev et all 



pie but powerfuU argument that the vast majority of stars 
must form in binaries: The lack of a significant single-star 
population in dynamically not evolved star forming regions 
means that stars cannot form in higher-order multiple sys- 
tems. These would decay on a system crossing-time (10^ yr, 
iGoodwin fc Kroupa|[2005l ). 

Af-body computations of initially binary-dominated 
star clusters have shown that a rising period BDF can be 
turned into a bell-shaped one by gravitational interactions 
among the clus ter members wi thin a few Myr (termed stimu- 
lated evolution, lKroupa|[l995al lbl. Marks, Kroupa fc Oh 2011, 
Oh et al., in prep.). Here, an analytical treatment for the 
change of orbital-parameter BDFs in star clusters, which ini- 
tially obeyed an invariant rising period BDF is used (Marks 
et al. 2011, hereafter Paper I) to efficiently calculate the 
stellar single and binary content in individual star clus- 
ters. Then the populations coming from all star clusters of 
a galaxy's freshly formed star cluster system are summed 
which yields galactic field stellar populations once the star 
clusters have dissolved (Dynamical Population Synthesis). 

We note that the method developed here to calcu- 
late the integrated galactic-field binary distribution func- 
tion (IGBDF, eq. [15] below) underlies similar concepts as 
the theory of the integrated galactic (stellar) initial mass 
function (IGIMF), which sums up the IMFs in all discrete 
star formation events showing that galaxy-wide IMFs are 
steeper at th e high-mass end than the invariant IMF in 
star clusters l|Kroupa fc Weidrieil l2003l : IWeidner fc Kroupal 
l2005h . This theory has proven extraordinarily successfuU 
in descri bing and pred ictin g observational properties of 
galax i es (iKoppen et al.ll2 007 : Pflamm-Altc nburg fc Kroupal 



20081: iPfiamm-Altenburg et al. .2009: Recchi et al 



Pfiamm-Altenburg fc KroupalbOOgTlCalura et ahlboid ) 



200' 



i 



A system refers either to a single star or to a binary. 



In Sec.[5]the model to calculate integrated galactic field 
binary distribution functions (IGBDFs) is devised and in 
Sec. [3] the results are presented and compared to observa- 
tions. Finally, Sec. 2] discusses and shows model predictions 
and Sec. [5] summarizes the main points of this investigation. 



2 MODEL 

In order to integrate over all stellar populations in discrete 
star formation events (Sec. I2.2ll . i.e. in embedded star clus- 
ters, it is necessary to first understand the evolution of bi- 
nary populations in them (Sec. I2.1|l . Note that the term dis- 
crete star formation event, star cluster and embedded star 
cluster are here used synonymously. 

In the following reference will be made to birth or pre- 
mam sequence and initial BDFs. The birth distributions 
(|Kroupal 1 1995al ) . including random-pairing of component 
masses for late type stars, describe the properties of bi- 
nary populations after they were born. But binaries are still 
embedded in their circum-stellar material, i.e. the compo- 
nents have not yet reached the ma in sequence stage. The ini- 
tial distributions (jKroupal Il995lf ) describe the correspond- 
ing statistical properties of a young binary population after 
birth binaries have undergone a phase of re-distribution of 
energy and angular momentum within their circumstellar 
material, called pre-mam sequence eigenevolution, acting on 
a time-scale of < 10^ yr. This mechanism introduces corre- 
lations between the orbital-parameters of short-period bina- 
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ries (such as between period and eccentricitv. lKroupall995bl ) 
as seen in observations. For a summary of these processes 
the reader is referred to Paper I. 



2.1 Binary distributions in star clusters 

In order to construct the field population by adding up the 
single stars and binaries in individual clusters, the evolution 
of binary populations has to be understood in terms of the 
initial properties of their host cluster (Meci,rh). 

In Paper I an efficient method is provided to analytically 
describe the first 5 Myr of the evolution of orbital-parameter 
BDFs in A^'-body computations of star clusters. The com- 
putations start with 100 per cent binaries distrib uted ac- 
cording to an initially rising period BDF derived in iKroupal 
l|l995bl ). which is consistent with constraints for pre-main 
sequence and Class I protostellar binary populations (see 
Paper I ). Their method qu antifies a stellar dynamical op- 
erator (|Kroupall200^ . 120081 ') . r^^yn "''''(*). in dependence of 
the initial cluster mass density, pcci = 3Meci/87rriJ, where 
Meci is the total mass in stars that formed in the embed- 
ded cluster and r^ is its initial half-mass radius. This op- 
erator transforms an initial {t = 0), perhaps invariant (see 
Paper I), orbital-parameter BDF, X>in, into an evolved one, 
^Af^^i.ri, j-^-j^ after some time t of stimulated evolution. 



«ccl,'-h 



(1) 



In the formulation of Paper I f^^y'n' ^^^^ particular on 



the initial BDF for binding-energies, <&iogio £5. in, 



(2) 



but extraction of other BDFs (period, semi-major axis, 
mass-ratio, eccentricity) is also possible with their model, 
given the interrelation between the orbital-parameters via 
Kepler's laws (Sec. 12.41 below). 

Let A'^cms = Ns + Ni, be the number of systems, i.e. 
the sum of all single stars and binaries with primary- 
star mass near mi in a stellar population. Then a BDF, 
^^eci.'-h^j^^-) _ logjjj iJ, a, e, g, . . .), is defined as the dis- 
tribution of binary- fractions, — Nt/Ncms, as a function of 



cE.f-''-»(mi) = ^ 
ax 



1 dNt{x) 



iVc, 



dx 



(3) 



where dNb{x) is the number of binaries in the interval 
[x,x + dx]. The total binary fraction equals the area below 
the BDF, 



fb{mi) = 



•(mi) dx 



(4) 



2.2 Integrated orbital-parameter distributions 

The galactic field's binary population is the sum over the 
populations in all star clusters that ever formed in a galaxy, 
having evolved for at least a time-span tfreeze after which the 
binary orbital-parameter properties become frozen-in, 



'■«col. 



dyn 



''■"(tfroezc) An Cocl(Mecl) rfM^d , (5) 



where ^eci(Afeci) is the embedded cluster mass function 
(ECMF). For simplicity it is assumed that all star clus- 
ters have formed with comparable half- mass radii, rh- The 
integration ranges from a minimum-cluster-mass, Meci,min 
to some maximum-cluster-mass, Med, max = Mcci,max(SFR). 
The maximum-cluster-mass that can form in a galaxy which 
has a given star for mation rate (SFR) is determined by 
(|Weidner et al.ll2004l '). 



M.d,max(SFR)^g^^^3 / SFR 



0.75 



(6) 



Star clusters are distributed according to a power-law 
ECMF with index /3, 



Cocl(Mc, 



X M, 



-/3 



(7) 



which is normalized such that the sum of the masses of all 
clusters equals the total mass of the freshly formed star clus- 
ter system, 

Msgs = / ' Med ecd(Med) dM^d . (8) 

The total mass, Msgs, needed to find the normalization con- 
stant fccd is determined from the SFR and the formation 
time-scale, 5t, of the star cluster system, 

Msgs = SFR x 5t . (9) 

IWeidner et"al] (120041 ) found that about every 3t = 10 Myr 
of ongoing star formation an ECMF is fully populated. 

The binary population will enter a galactic field with 
characteristics set at time tfrcozo- This is reached after suf- 
ficient time for cluster internal stimulated evolution when 
only hard binaries are left (Paper I) , or if a cluster suddenly 
expands rapidly, e.g. as a result of residual-gas expulsion, 
which inhibits further stimulated evolution. This state is al- 
ready reached after a few Myr or even in less than 1 Myr for 
de nse configuration s llKroupa 1995a: D uchcnc ct al. 19991: 
iFregeau et ahl l2009l : iParker et all [2OO9I . Paper I). For our 
purpose, we choose the freeze-in time to coincide with the 
time of the occurence of the first supernovae, tfrcozc ~ 3 
Myr. If the cluster is still embedded, supernovae are ex- 
pected to drive out the residual-gas of the embedded cluster 
rapidly, leading to cluster expansion and destruction of the 
maj ority of clusters in a star cluster system («90 per cent, 
e.g. lLada fc Ladall2003l ). The exact time star clusters are al- 
lowed to evolve their population is not that important since 
the time-scale on which binary dissolution occurs is short. 
In particular, the difference in the binary fraction between 
3 and 5 Myr in the A''-body models used in Paper I is of the 
order of a few per cent only. 



2.3 Normalization 

Since the stellar-dynamical operator (eq. 01 transforms be- 
tween BDFs independently of the number of systems in a 
star cluster, one has to take care to preserve the definition 
and normalization for individual clusters (eqs. [3]and|3| also 
for the integrated population. Consider therefore as an ex- 
ample two clusters consisting of A'cms = 10 and 100 systems, 
with ft = 40 and 60 per cent, respectively. Evaluating eq. ([SJ 
for the two clusters would result in ft — 50 per cent for the 
combined (or integrated) population, but the true resulting 
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population has ft, — 64/110 = 58 per cent. Thus, the num- 
ber of systems making up the population has to be taken 
into account. 

Define a BDF for a galaxy analogous to eq. (|3]), 



ax 



dx 



(10) 



where dA^f ^(a:) and N^Is 

are the number of binaries in the 
interval of size dx and the number of cm-systems in a whole 
galaxy, respectively. The goal is thus to calculate the num- 
ber of binaries per a;-interval and number of systems (sin- 
gle-|-binary) in a galaxy from the respective numbers in star 
clusters separately. 

For an initial cluster mass, M^d, the total number of 
freshly hatched stars in that particular cluster is calculated 
from 

Mecl 



iV.(^ccl) = 



(11) 



where m ~ OAM^) is th e average mass of the canonical stel- 
lar IMF (|Kroupall200ll ). The total number of binaries is re- 
lated to A''.(Meci) via 



Nt(Mccl) = Af,(A/ecl) - iVcms(A/ocl) 



(12) 



Inserting this in — Nt/Ncms and rearranging gives the 
number of cm-systems in that particular cluster, 



A*'cms(Mccl) = 



(13) 



Since $ 



is known (eq.[2} and therefore (eq. |4]), 
A'cms and Nt can be calculated. The number of binaries per 
interval dx becomes. 



^^^^^^=iV_(A/e,)x<i>f-'-^. 
ax 

Therefore, eq. (|10p becomes. 



<I>?^(mi) = 



dx 



Cecl(Mecl) dMecl 



(14) 



.(15) 



Eq. (jlSp is referred to as the integrated galactic-field binary 
distribution function (IGBDF). It can be calculated for sin- 
gle stars and binaries with primary-mass in an interval Am 
around mi, e.g. 0.8 — l.O4Af0, or for all late-type stellar 
systems mi/Af© £ [0.08,2] (Tab.©. 



2.4 P-, e-, q- and a-IGBDF 

In order to extract IGBDFs for different orbital-parameters 
(period, P, eccentricity, e, mass-ratio, q, and semi-major 
axis, a) from the known energy IGBDF, the procedure is 
basically as described in sec. 4.2 of Paper I for orbital- 
parameter distributions in single clusters. The idea is to 
compile a large library of binaries (A'lib = lO'' binaries for 
the present purpos e) whose proper ties are selected accord- 
ing to the recipe in iKroupal l|l995bl . see Paper I). Here, the 
library consists only of binaries with primary masses up to 
2Mq in order to resemble a galactic field and it contains 
the values for mi, m2, Eb, P, e, q and a. The total number of 
initial binaries in the integrated population, calculated ac- 
cording to Sec. 12.31 is scaled so as to match the library size 
of Niih binaries. Following this, the final distributions are 



Table 1. Adopted mass-ranges for single-stars or primary- 
components of binaries with different spectral- type (SpT). 



SpT 


F 


G 


K 


M 




1.04-1.4 


0.8-1.04 


0.45-0.8 


0.08-0.45 



constructed by removing an appropriate amount of binaries 
with a given binding-energy from the library according to 
the number-ratio of binaries in the resulting and initial en- 
ergy IGBDFs. The remaining binaries are used to construct 
the P~,e~,q— and a— IGBDF. In order to extract sub- 
distributions, such as for a special spectral-type or period- 
range, only those remaining binaries are used which fuUfiU 
the additional criteria. In the forthcoming sections, mass- 
ranges for single-stars or primary-component masses for 
binary-stars, respectively, are adopted, as shown in Tab.[T] 



3 RESULTS 

According to eq. (|15|l . the properties of an IGBDF depend 
on the minimum embedded cluster mass and the SFR (i.e. 
the maximum cluster mass, eq. [6| of the considered galaxy, 
the index, /3, of the ECMF (eq. [7]), and the average half- 
mass radius, rh, of the star clusters. The influence of these 
parameters on the integrated binary properties by means of 
the global binary fraction will be investigated. A discussion 
of the Galactic field binary population of the Milky- Way 
(MW) and what can be learned about the initial star cluster 
system of the MW from the observed binary population ends 
this section. 



3.1 The importance of low-mass clusters for the 
Galactic field binary population 

From which type of cluster do most Galactic field binaries 
originate? Dissolving low-mass embedded clusters will each 
contribute only a small number of systems to a Galactic field 
population but a large fraction of binaries will be among 
them due to inefficient stimulated evolution (Paper I). For 
high-mass clusters the situation is the other way round. 
However, the ECMF (eq. (Tjl is steep, i.e. there are many 
more low-mass than high-mass clusters. Thus, this question 
is not simple to answer qualitatively. 

Therefore the number of systems, binaries and single 
stars which are added to the Galactic field by all clusters of 
a given mass Med is calculated by evaluating 



Af„ns/b/s(Mccl) Cccl(Afccl) . 



(16) 



dMsd 

The result is depicted in the left panel of Fig. [T]for the MW 
star cluster system with j3 — 2.0, rh — 0.2 pc, Afcci.min = 
5Mq and SFR^SA/q yr"^ (see Sec. [33] below) . It shows 
that the steepness of the ECMF dominates such that most 
systems in the Galactic field stem from low-mass clusters. 
For the same star cluster system parameters, the right panel 
of Fig. [T] plots the normalized cumulative number. 



cms/fe/s 



dMcci 



dMe, 



(17) 



where A^c^^fs/b/^ 



is the total number of systems, binaries or 
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10" 10^ 10^ 10^ 10"* 10^ 10^ 10" 10^ 10^ 10^ lo'* 10^ 10'' 

cluster mass M^^^ I Mg^r, cluster mass Mg^i / Mgu^ 



Figure 1. Left panel: Number of all singles (dotted), binaries (dash-dotted) and all systems (single+binary, solid) per cluster mass in 
the Galactic field (eq. ll6ll , respectively, that form within one star cluster system formation time-scale <5t (eq.|9]l. Lines are drawn for the 
MW star cluster system (/9 = 2.0, rjj = 0.2 pc, Mgei min = 5M0 and SFR=3M0 yr"-*^). Low-mass clusters contribute most systems to 
the Galactic field population owing to the steepness of the ECMF. The intersection between the single-stars and binary-stars line is the 
cluster mass at which = 50 per cent for the used parameters. Right panel: Cumulative number of singles, binaries and all systems 
as a function of Mgd normalized to the respective total number (eq. I17II for the same star cluster system parameters as in the left panel. 
Low-mass clusters (Mod li, SOOMq) are the dominant binary contributors while high- mass clusters (Af^ci ^ W^Mq) donate most single 
stars. 



singles in the Galactic field, respectively, i.e. eq. (I16p in- 
tegrated from Mod, min to Med, max- It demonstrates tiiat 
about 50 per cent of all systems in tiie Galactic field formed 
in clusters with initial masses Mod < 2 X 10^ M0, while 
50 per cent of all Galactic field binaries originate from clus- 
ters witli Afoci % 300Mq only. For single stars the situation 
is inverted. Roughly 50 per cent of all single stars come from 
star clusters witii masses Mod ^ 10'* Mq. 



3.2 Parameter study &: degeneracy 

Fig. [2] depicts the global binary fraction (color-coded), i.e. 
including all late-type binaries (mi ^ 2Mq), of a galax- 
ies' stellar population in dependence of the four parameters 
(Meci.min, SFR, /3, Th) in eq. ([TSjl . We learn the following from 
tliese diagrams: If all other parameters are fixed... 

... the higher the SFR, the smaller the field-binary popu- 
lation. A higher SFR yields a larger Mod (eq.[6]). A higher- 
mass cluster is generally denser than a lower-mass cluster 
and will therefore contribute a population with a lower bi- 
nary fraction (Paper I). Therefore drops with increasing 
SFR. 

... the larger rh, the larger the galactic field binary frac- 
tion. Larger typical half-mass radii imply lower densities, 
resulting in dynamically less evolved binary populations be- 
fore dissolution (higher ft, Paper I). Therefore ft increases 
with increasing rh. 

... the steeper the ECMF (larger index j3), the larger the 
field-binary population. A steeper ECMF amounts to a lesser 
proportion of high-mass, i.e. higher-density clusters. There- 



fore on average populations with larger binary-fractions en- 
ter the field and increases. 

... the larger the mimmum cluster mass, the fewer bina- 
ries exist in the field. Increasing Afod.min means cutting out 
of low-density clusters which would contribute populations 
with large binary proportions. The remaining higher-mass 
clusters shed dynamically more evolved populations into a 
galactic field. 

Fig. |2] distincly shows that the parameter space is degen- 
erate. E.g., the increase in ft due to a steeper ECMF can 
be compensated for by introducing smaller cluster sizes (left 
panel) or assuming a larger A/oci,min (right panel). Similarly, 
a larger typical rh can be counteracted with a higher SFR 
(middle panel), and so on. Therefore constraints on at least 
some of the four parameters from independent sources are 
needed to reduce the allowed solution space, when compar- 
ing model predictions with observed distributions. Those are 
available for the MW. 



3.3 The Milky-Way 

In order to calculate a synthetic stellar population for the 
MW's Galactic field, we need to estimate the parameters 
entering the IGBDF for our Galaxy. 

The current global SFR of the MW using different 
met hods is d e termin e d to lie betw e en » 1 and 5M0 yr~^ 
(Sm ith et all Il978l: iDiehl et al.1 I2OO6I: iMisiriotis et all 
200fi ICalzetti etaT 



20091 : iMurrav fc RahmanI I2OI0I : 



iRobitaille fc Whitney! I2OIOD . These values are compatible 



with the assumption that the total ste llar mass in the disk 
and bulge of the MW (5 x 1O^''M0, iBinnev fc Tremaind 
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Figure 2. Study of the influence of the parameters in eq. {TS} determining the integrated stellar population. In each panel the color- 
coding shows the global binary- fraction, /(,, i.e. including all late-type binaries (mi ^ 2Mq), according to the bar on the right of each 
panel (in per cent). The solid overlaid lines indicate curves of constant binary fraction. The big solid white dot in each panel marks the 
position of the MW (Sec. |331| . The binary population increases with decreasing SFR, with increasing cluster half-mass radii, rji, with 
steepening of the ECMF (increasing index /3) and with decreasing minimum cluster mass, M^^-j ^j^. 



l2008h has assembled continously during the last 13.7 Gyr 
with a SFR of S.GM© yr"^ Assuming SFR^SAf© yr~^ 
yields Meci.max = 1-9 X 10^ A/0 (eq. [6| for the most 
massive cluster having formed in the M W disk, compa- 
rable to the most massive open cluster in iPiskunov et al.l 
120071 ') 's sample of 236 open clusters within 1 kpc from 
the Sun (1.1 x 10^ Mq for Sco 0B5). Taking the min- 
imum clust er mass similar to that of a Taurus- Auriga 
like group (|Kroupa fc Bouvieij l200l A/cci,min = 5Mq) 
and an ECMF index /3 = 2 similar to the observed 
slope s for Galactic embedded clusters (iLada fc Ladal 



20031: Id e la Fue nte Marcos fc de la Fuente Marcoj l2004l : 



Gieles et al. 2006), we can evaluate the properties of the 



Galactic field binary population by numerically integrating 
eq. (|15[) for different embedded cluster half-mass radii. 



0.3 



0.25 



LU 

? 0.2 

o 



LU 

o 



C3) 
O 



0.15 



0.1 



0.05 



r,, / pc'= 8^ 
ru / pc = 0.6 

i-h/pc = 0.2 
rj, / pc = 0. 1 



3.3.1 Energy distribution 

Fig. shows the energy IGBDF for the above parameters, 
for different cluster half-mass radii and for primaries of all 
masses. It is found that the binary-fraction in the Galactic 
field (the area under the distributions) is smaller if the typ- 
ical cluster radii are smaller. Smaller average radii lead to 
higher initial densities in star clusters which results in more 
efficient binary dissolution and thus a lower binary-fraction 
(Sec. 13. 2[) . For rh = 0.1 pc the overall binary- fraction is 
fb = 0.34 and for rh = 0.8 pc, ft = 0.71. 

3.3.2 Period distribution 

For the adopted initial conditions in star clusters (Paper I) 
and the above values for the SFR, /3 and Afeci.min in the 
MW, the period IGBDFs in Fig. g] for different values of rh 
demonstrate that the initially rising period distribution in 
star clusters translates into a bell-shaped form in the Galac- 
tic field, at least if the star clusters from which the Galactic 
field population originates, were rather compact. Compari- 
son with the corrected observed Galactic field period BDF 
for G-dwarfs (DM91) suggests that clusters formed quite 
compact (left panel). The G-dwarf field binary fraction of 
~ 57 per cent is best reproduced with a typical initial cluster 
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Figure 3. Global energy IGBDFs for different typical initial clus- 
ter half-mass radii in the MW. The smaller the clusters are, the 
larger are their densities and the more efficient is the binary 
destruction initially (Sec. 13.21 1. Therefore the binary fraction is 
smaller the lower the initial value of rh is, ranging from = 0.34 
for rh = 0.1 pc to /j, = 0.71 for rh = 0.8 pc. The area under each 
distribution equals the total model binary-fraction in the Galactic 
field. 

size of rh = 0.1 p c, wh ere fb{G) = 0.58. A similar study by 
iRagh avan et all (|2010l ) investigating the multiplicity prop- 
erties of 454 solar-type stars selected from the Hipparcos 
catalogue is in agreement with the DM91 data and thus 
wit h the compact formation of star clusters. 

iFischer fc Marcvl (|l992l ) studied the M-dwarf binary 
population within 20 pc from the Sun, showing that the pe- 
riod distribution of M-dwarfs can also be described by a log- 
normal distribution. Comparison with the IGBDF models 
(Fig. |4l right panel) favors their formation in slightly more 
extended clusters. For rh = 0.3 pc the binary-frequency 
of 41 per cent compares well with the observational result 
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Figure 4. Period IGBDFs for the MW model (Sec l3.3l l with different typical initial clus ter radii for G-dwarf (l eft) and M-dwarf (right) 
binaries (histograms). Comparison with the G-dwarf period BDF by DM91 (circles) and lRaghavan et al. I | |201C| . triangles) suggests that 
clusters should have formed very compact. The period I GBDF with a ha, l f-maas radius of = 0-1 pc matches t he observed binary- 
frequency best. The observed M-dwarf binary- frequency jFischer fc Marcvl lll992l . robust) and lBergfors et al.l 1 I2OIOI) ) favours formation 
in = 0.3 pc sized clusters. This difference in rji between the G-dwarf and M-dwarf solution is discussed in Sec. U] The line-types are 
as in Fig. |3] 



(/i, = 0.42). Recentlv. lSergfors et all l|2010h surveyed 124 Ad- 
type stars for binary separations a < 200AU for their multi- 
plicity properties. The observed semi-major axis distribution 
has been translated into periods using Kepler's laws and are 
incorporated in Fig. [4] (right pane l ). Th eir data are com- 
patible with the iFischer fc Marcvl (|l992l ) data. For a sub- 
sample of the model binary population with a < 200 AU, as 
in the observations, the rh = 0.2 pc-model binary-frequency 
is /i, = 0.31 and ag rees with the o b serve d binary-fraction 
of ~ 32 per cent. iBergfors et al] l|2010h note, however, 
that there might be some overabundance of systems with 
P < 20 d, corresponding to a < 80 AU for an average system 
mass of « O.5M0, due to the sample selection and that cor- 
rections might be necessary for non-physical (optical) pairs. 



3.3.3 Mass-ratio distribution 

The result that star formation in rather compact structures 
is favoured is independently confirmed by considering the 
distribution of mass-ratios. Fig. [5] depicts the mass ratio 
IGBDFs for the complete field population (left, IGBDF for 
all binaries with primary masses mi 5C 2i\f0), for G-dwarf 
binaries only (middle) and for M-dwarfs (right). It is appar- 
ent that the complete distribution is fiat between roughly 
q = 0.2 and 0.9 while the G-type binary sub-distribution is 
decreasing and the M-dwarf IGBDF is increasing with in- 
creasing q. We emphasize, that all three mass ratio IGBDFs 
result from initially sampling the two birth components of a 
binary random ly from the same underlying stellar IMF (e.g. 
lKroupa|[l995al '). The decreasing trend is also evident for F- 
and K-type primaries (Fig.[S]below) and is seen in the obser- 
vations of G-dwarfs by DM91 (middle panel of Fig. [5]). The 
observational data compare very well with the compact for- 



mation models, as above for the period IGBDFs. Addition 
of the mass ratio IGBDFs for the different spectral types 
(Tab. [T] and Fig. [S] below) results in the flat distribution for 
the complete mass ratio IGBDF (left panel in Fig. [5]). The 
binary fraction of 35 per cent in the observations o f 106 G - to 
M-type systems in the analysis by iReid fc Gizij l|l997l ) are 
best compatible with the rh = 0.1 pc model where fh = 0.34. 

In contrast to the declining q-distribut ion for G-dwarfs 
found by DM91, note that the results of iRagh avan et al.l 
(|2010t ) suggest a mass-ratio distribution for binaries with a 
solar-type primary which is fiat between q « 0.2 and 0.9, 
while fi nding a similar per i od di stribu tion (see discussion in 
Sec. BJ. iFischer fc Marcvl (|l992l ) and iBergfors et al] (|2010l ) 
extracted also a mass ratio distribution for their respective 
samples of M-dwarfs whose rising shapes are well reproduced 
by the models (right panel). As expected from the period 
IGBDFs (Fig. |4|, the same rh = 0.3 and 0.2 pc models, re- 
spectively, are consistent with the observations, being again 
larger than for G-dwarfs. 



3.3.4 Eccentricity distribution 

The eccentricity IGBDF for G-dwarf binaries is in agree- 
ment with the observational data (Fig. |6]|. The distribu- 
tion of eccentricities is bell-shaped for logj^Q Pcirc ~ 1-06 < 
logjQ P < 3, where Pdrc is the circularization period identi- 
fied by DM91 below which the orbits are circular (e f» 0). 
Circula risation occurs through pre-main sequence eigenevo- 
lution ijKroupal [T995bl . Paper I). The eccentricity IGBDF 
follows the thermal distribution for logjQ P > 3 because it 
is invariant to stimulated evolution. 
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Figure 5. Mass-ratio IGBDFs for all latc-type binaries (mi ^ 2Mq, left panel), G-dwarf (middle panel) and M-dwarf binaries (right 
panel), for the MW model (Sec. [33)l with different half- mass radii (histograms). The line-types are as in Fig. [3] While the distribution 
between q = 0.2 and 0.9 is flat for the complete distribution it declines for the G-dwarfs and increases for M-type binaries with increasing q. 
The peak at mass ratio's in the model close to unity (an effect of pre-main sequence eigenevolution) is also evident in the o bservational 
data for all primary-masses combined jReid fc Gi2ij|l997l . left) and M-dwarf binaries jPischer fc Marcvlll992l :[ Bergfors et al l20ia . right), 
while it is less pronounced, if at all, for the G-dwarf observations by DM91 (middle). Note that agreement is obtained for the similar rji 
as for the period IGBDFs (Fig.lDl. 
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Figure 6. Eccentricity IGBDFs for G-dwarf binaries in comparison with the observations by DM91 in the indicated period ranges, 
^circ = 11-6 d is the circularization period derived from the observational data. Note that in both panels the data is (exceptionally) 
normalized to the total number of binaries instead of systems, in order to be able to compare to the observations. This is also why IGBDF 
models with different (histograms, different line-types) can hardly be distinguished. Left panel: The eccentricity BDF is bell-shaped 
for orbital periods below 10^ d due to pre-main sequence eigenevolution as in the observational data. Right panel: For P ^ 10^ d the 
e-distribution follows the thermal distribution (/i,(e) = 2e, solid line) for the IGBDF model as well as in the observations. The thermal 
eccentricity distribution is invariant of stimulated evolution. 



3.3.5 Single and binary population in dependence of 
spectral type 

The fraction of binaries in the field is a function of the spec- 
tral type (i.e. mass) of the primary. The later the primary 
spectral type the lower is the resulting binary fraction in 
the field (Fig. [7|). Especially M-dwarfs have a lower binary- 
fraction than binaries with a F-, G- or K-type primary. The 
reason for this is twofold. 

First, due to the shape of the stellar IMF M-type stars 
are most numerous so that upon random pairing of the bi- 



nary components at birth and after eigenevolution (Sec. [2| 
about « 90 per cent of all initial binaries carry a com- 
panion of spectral type M. Thus, everytime a binary dis- 
solves, in 9 out of 10 cases at least one M-dwarf will end 
as a single star (two if a M-dwarf binary dissolves). Each 
M-dwarf contributes to A'"cms(M), and therefore /^(M) — 
Nb{M)/N,ms{M) shrinks. 

Secondly, binaries with a lower binding-energy, Ei,, are 
more prone to dissolution (Paper I) and the binding-energy 
is proportional to the mass of the primary, Et oc mi. Since, 
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to 55 — 66 per cent and is in excellent agreement with the 
estimate bv lLadal (|2006l ). that about 2/3 of all primary stars 
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Figure 7. Comparison of the binary fractions among systems 
of one spectral type (SpT) for IGBDF models with different 
(filled and open circles for an ECMF slope /3 = 2 and 2.4, re- 
spectively, in each column from top to bottom: 0.8,0.6,0.4,0.2 
and 0.1 pc). The initial binary fraction for all spectral types is 
/(, = 1 (dotted horizontal line). M-dwarf binaries have the low- 
est binary-fraction when they enter the Galactic field. The flatter 
the ECMF (smaller 0), the lower is the binary-fraction (Sec. l3.2ll . 
Observed binary fractions agree with the /3 = 2, = 0.1 — 0.2 pc 
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by construction, all binaries intially follow the same initial 
period distribution (jKroupal [l995bl . Fig. [4|, the M-dwarf 
(low mi) energy IGBDF is shifted to slightly lower energies 
compared to F-, G- and K-type binaries. Thus, it is gener- 
ally easier to dissolve M-dwarf binaries than binaries with a 
primary of an earlier spectral type. Furthermore, low-mass 
binaries are more frequent, again owing to the shape of the 
stellar IMF. In the model, ~ 57 per cent of all binaries have a 
M-dwarf primary initially. Dynamical encounters including 
M-dwarfs will therefore occur often in the clusters. In turn, 
each disruption of a M-dwarf binary will reduce Ni,{M) by 
one, increase A'^cms(M) by one and, in turn, reduce /i,(M). 

All observational data in Fig.[7]again compare best with 
star formation in compact star clusters (rh = 0.1—0.2 pc and 
/3 = 2). Although all stars are locked up in binaries initially, 
more than half of all systems in the Galactic field end up 
as single stars due to stimulated evolution in star clusters 
before they dissolve (Fig. [2]). M-dwarfs in the IGBDF model 
constitute about ~ 80 per cent of the single star population 
in the Galactic field. Only ~ 13, 3 and 2 per cent of all single 
stars in the IGBDF model have spectral type K, G and F, 
respectively (the remaining 2 per cent are of spectral type 
A, for a Galactic field population which consists of systems 
with mi ^ 2 Mq, Sec. [23]). 

For rh ~ 0.1 — 0.2 pc and P — 2 the single star fraction, 
/s(SpT) = 1 - /i,(SpT), in the model becomes f,{M) ^ 
0.75 - 0.85, fs{K) f» 0.36 - 0.48, fs{G) ^ 0.31 - 0.42 and 
/s(F) ~ 0.28 — 0.39. The total single star fraction amounts 



are single. 



4 DISCUSSION & MODEL PREDICTIONS 

4.1 Formation in compact star clusters 

Comparison of the observational data with the model has 
suggested that MW star clusters typically formed quite com- 
pact (rh ~ 0.1 — 0.3 pc. Sec. 13. 3p . Such small radii compare 
with the observational lower end of the sizes of dense cores 
in giant molecular clouds. However, the sizes of these dense 
cores range up to 2 pc and the spatial extends of embedde d 
clusters are typically comparable (e.g. iLada fc Ladall2003l ). 
But there is evidence that the forming s tars within the em- 
bedded cluster sta r t dynamically co ld (|Walsh et al.l |2004| : 
iPeretto et alll2006l : iLada et al.l 120081 ') and that protostellar 
objects are more confined than more evolved yo ung stellar 
objects (|Teixeira et al.ll2006l : iMuench et al.ll2007l ). Thus, an 
embedded cluster will collapse to a smaller configuration. 
In that sense the here derived half-mass radii may be inter- 
preted as the sizes of clusters when they reach their peak 
stellar density and stimulated evolution is most efficient. 

Additionaly the IGBDF model implicitely assumed that 
a typical half- mass radius for all clusters exists (Sec. 12. 2p . 
If a mass-radius r elation for embedded c lusters better de- 
scribes reality (e.g. [Harris fc Pudrit j[l993 . for virialised gas 
cores in giant molecular clouds), the typical rh for MW clus- 
ters can be seen as an average value for all clusters of any 
mass and size. Note that if such a mass-radius relation ex- 
ists, upon averaging the inferred best rh will be closer to 
the true value for low-mass clusters, which are most im- 
portant for the Galactic field binary population (Sec. 13. l|) . 
Higher mass clusters will then occupy a somewhat larger ra- 
dius range. However, a trend of radius with luminosity (or 
mass) in young star clusters has been shown to b e shallow 
for the galaxy merger NGC 3256 (Zcv>f et al."l999'), for stel- 
lar clusters in 18 spiral galaxies (|Larsen 200 4) and for young 
clusters in M 51 jScheepmaker et al.ll2007l '). 

The typical rh is arrived at assuming P = 2 down to 
5Mo. The ECM F might however fl atten (/3 0) below 
f» 50 - IOOM0 (|Lada fc Ladal 120031 ') . implying that fewer 
low-mass clusters are present compared to the numbers used 
in the model. According to a computation with a broken 
power-law ECMF, where /3 = for Med lOOM© and 
P — 2 otherwise was chosen, typical half-mass radii nedd to 
be larger by « 0.1 pc in order to agree with the observations. 

A small difference in the solutions for Hi (0.1 vs. 0.3 pc) 
in the solutions for G- and M-dwarf binaries might be evi- 
dent (Fig. [4]). Since the estimated cluster size is also a mea- 
sure of how dense the region is in which the respective sub- 
population has formed, the possible difference might simply 
indicate their formation in different locations of the same 
cluster. A primordial mass-segregated cluster, where the G- 
dwarfs would be more centrally confined to a region of higher 
density while M-dwarfs form out to larger radii, would nat- 
urally account for the apparently different ranges of allowed 
cluster radii. 

The uncertainty in the inferred typical cluster size of 
0.2 pc does not appear to be very large, given the estimates 
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for rh by comparison of the model with independent observa- 
tions (Sec. I3.3|l . Even if a poss ible error in t he ECMF index, 
(3, of up to 0.5 is considered l|Larsenl ^2009) the uncertainty 
in rh is of order 0.1 — 0.3 pc only (Fig. (2] left panel). 

To calculate the composition of the Galactic field stellar 
population a (constant, average) global SFR was adopted. 
A declining SFR history might instead be better suitable to 
descr i be the evolut ion of the MW di sc dSoissier fc Prantzo^ 
1 19991 : iNaab fc Ostriker 2006 ; Schonrich fc Binnevll2009l ). If 
early on the SFR has been higher than average, clusters 
more massive than allowed for the adopted SFR (eq. [S]) 
would have been able to form, which would have contributed 
a larger number of single stars. Later, when the SFR sank 
below the average SFR the highest mass clusters wouldn't 
be able to form any more and more binaries would en- 
ter the field originating in the low-mass clusters. These ef- 
fects might eventually compensate each other, but will de- 
pend on the actual history, which can't be tested without 
a modification of the used code. Also, the global history 
might be non-representati ve for the solar neighbourhood 
l|Boissier fc PrantzosI IT999I ) . However, even if the SFR has 
been much higher in the past only and settled to the present 
value, the inferred typical rh wouldn't change strongly. Re- 
ducing or enhancing the SFR by up to two orders of mag- 
nitude and at the same time retaining the observed binary- 
fraction requires, respectively, a smaller by ~ 0.1 pc or 
larger by « 0.2 — 0.3 pc only (as evident from Fig. (2] middle 
panel) . 

In this sense, the typical cluster size seems to be well 
constrained by the IGBDF model. 



4.2 MW orbital-parameter BDFs 

Since the IGBDF models for the MW (Sec. [33} with typical 
star cluster sizes of rh « 0.2 pc agree well with independent 
observational data, these models are used to predict the pe- 
riod, semi-major axis, mass-ratio and eccentricity IGBDFs 
for the solar neighbourhood. 

Fig. [8] depicts the resulting distributions for binaries 
of different spectral type and the combined distributions 
(mi ^ 2Mq; no additional cuts, e.g. in period, are applied). 
For the period, mass-ratio and eccentricity IGBDFs the dis- 
tributions for F-, G-, and K-type binaries are very similar 
and therefore probably hard to distinguish by observations. 
However, the mass-ratio IGBDF might hold the ability to 
test the model prediction for different spectral types at low 
(g « 0.1 — 0.2) and high {q « 0.9 — 1) mass-ratios since dif- 
ferences are more pronounced there. In particular, there is a 
larger gap between the G- and K-binaries around q « 0.15 
and between the F- and later type binaries at g f» 0.05 which 
might be visible in observations, too. 

The M-dwarf binary IGBDFs for all quantities are dis- 
tinct from the corresponding distributions for earlier spec- 
tral types since breaking-up of binaries having at least one 
M-type component happens frequently before their birth 
clusters dissolve f Sec. 13. 375]) . Thus, by the time the M-dwarfs 
emerge from the clusters their major ity are single stars con- 
sistent with observations (lLadall200^ . Sec. 13. 3. 5)) . despite be- 
ing born as binaries. Differences in observed distributions for 
M-dwarf binary populations and those of later types should 
be apparent and are thus suited to test the IGBDF pre- 



dictions. The combined IGBDFs lie typically between the 
M-dwarf distribution and the earlier types. 

Note that it is not possible to show a common initial 
mass ratio distribution since it depends, in contrast to the 
other IGBDFs, on the c onsidered spectral type (see also 
iKouwenhoven et al.ll2009l V This is a result of the random 
selection of birth binary component masses and the upper 
limit for the mass of the primary star when a spectral type- 
limited sample is investigateqj. While at least the shape 
of the initial mass ratio distributions is similar for F-, G- 
and K-binaries it is completeley different for M-dwarfs (com- 
pare, e.g., the middle with the right panel in Fig. [5] for the 
initial distributions of G- and M-dwarfs, respectively). The 
IGBDFs decline with increasing mass ratio for F-, G- and 
K-type binaries, while the mass ratio IGBDF for M-dwarfs 
is increasing with increasing q. We explicitly note that the 
same observational data for M-dwarfs, for which typically a 
flat g-distribution is inferred, is consistent with the increas- 
ing trend in the models. The mass ratio distribution is fiat 
only if all primary-masses are combined to construct a mass 
ratio IGBDF. 

As opposed to DM91, iRaghavan et al.l l|2010l ') find a 
mass ratio distribution for binaries with a solar-type primary 
that is fiat between q ~ 0.2 and 0.9, while finding a simi- 
lar period distribution as DM91. This can not be expected 
within the framework of the IGB DF model. Since the ran ge 
of considered primary masses in iRagh avan et al.l (|2010l )'s 
and DM91's study are comparable, their mass ratio distribu- 
tions should be in agreement if the period distributions are, 
and vice versa. A flat distribution from the models is only 
obtained if the complete late-type binary population is con- 
sidered (Fig. [5] left panel, and Fig. [8]). Although this appar- 
ent inconsistency might be due to model assumptions (see 
Sec. l4.4|) . the presented results are otherwise very successfuU 
in describing independent observational data for the MW, 
adjusting only one free parameter (the typical half-mass ra- 
dius rh, Sec. I3.3|l . It is noted that DM91 monitored radial 
velocities for their sample over a p eriod of 13 years find - 
ing accurate orbital solutions, while IRaghavan et al.l (|2010l ) 
compile data from various sources covering many tech- 
niques. However, in order to understand these issues better 
it will be neces s ary to study t he consistency betw een the 
IRaghavan et all l|2010l ) and the iReid fc GizTdi |l993) data - 
how can both, the (essentially) G-dwarf and the all-primary 
combined mass ratio distributions be flat at the same time? 



4.3 Dependence on galaxy morphology 

The binary properties of galaxies within the IGBDF model 
are strongly dependent on the SFR fSec. 13.2"]) . Since SFRs 
in galaxies are observed to cover a large range from « 
Mq yr"^ for giant eUipticals to > 1 000 Mp yr~^ for uhra 
luminous infrared galaxies (ULIRGS. lGrebel l l201ll ), binary- 
frequencies and binary-properties are expected to vary be- 
tween galaxies of different morphology. This might have cos- 
mological implications, such as for the SN type la rates in 
galaxies. 



By definition the secondary mass can only be lower. 
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Figure 8. Predictions for the period (logjQ P), semi-major-axis (a), mass-ratio (g) and eccentricity (e) IGBDFs for different spectral 
types and the combined distributions in the solar neighbourhood for the MW models (Sec. 13.311 which best fit the observational data 
(rji 0.2 pc). The logjQ P-, a- and e-IGBDFs for F-, G- and K-binaries appear very similar, but can probably be distinguished at low 
and high g- values in the mass-ratio IGBDF. The peaks seen for q = 0.9 — 1 in the mass ratio distributions are due to pre-main sequence 
eigenevolution (Sec. O. All IGBDFs for M-dwarfs and the all-primary mass IGBDFs are distinct from those for F- to K-binaries. Note 
that a common initial distribution for the g-IGBDFs can not bo given since they are different for different spectral types (see discussion 
in Sec. 



4-3.1 Elliptical galaxies 

Ellipticals (Es) are nowadays more or less free of cold gas and 
are pressure- or random-stellar-motion-supported with low 
or no star formation activity. Even the suspected precursors 
of giant ellipticals, quasars at high redshift {z 6, tuniverse < 
1 Gyr), reveal supersolar metallicities indicating a starburst 
that quickly enriched the material with metals l|Fan et alj 
[2003). This suggests that Es had a large SFR initially until 
their gas reservoir was depleted. One of the highest-redshift 
quasars kno wn has a S FR of ~ IOOOMq yr~^ as derived 
for ULIRGS l|Fanll2006D . For SFR= 10^ Mq yr-\ from the 
middle panel of Fig. [2j a binary fraction of the order « 
30 — 40 per cent for a typical cluster size of rh — 0.2 pc can be 



inferred. It is noted that a non-shallow mass-radius relation 
(Sec. I4.1|l . which is not considered in the present models, 
might affect the results for such high SFR. If, additionally, 
during star bursts low-mass clusters are not able to form, 
i.e. Mcci,min IS larger, this would further lower the binary- 
fraction in Es (Fig.[2l right panel). Thus, if E galaxies formed 
rapidly they ought to have low binary fractions. 



4.3.2 Spiral galaxies 

In terms of the SFR, spiral galaxies are intermediate objects 
between Es and dwarf galaxies. SFRs in spirals like the MW 
lie between 0.1 and IOMq yr"^ (e.g. iLee et "ai1l2009l . [ioill '). 
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Figure 9. Predictions for tlie period (logj^Q P) and mass-ratio (g) IGBDFs in dependence of galaxy morphology. Model parameters are 
as those for the MW (Sec. 13.31 rjj = 0.2 pc) except that SFRs are used according to the galaxy type (Sec. 14.31 1. Elliptical (E) galaxies 
have the lowest binary-fraction (the area below the distribution) and the period IGBDF peaks at shorter periods than for spirals (S) 
and dwarf irregulars (dIrr). The shapes of the mass ratio IGBDFs for the different galaxies resemble each other. The pre- main sequence 
eigenevolution peak in the mass ratio IGBDF is visible, as in Fig. |8] 



For the same cluster size and a SFR of IMq yr~^, from 
Fig. [2] a global binary frequency of « 40 — 50 per cent is 
expected, similar to what is seen in the solar neighbourhood 
(Sec. [33}. 



4-3.3 Dwarf irregular galaxies 

Dwarf irregular (dIrr) galaxies have very low SFRs rang- 
-5 n/r_ ,,^—1 



ing down to 10"^ Mq yr"^ (|Lee et al.lT2009l . 1201 ll ). For 
rh ~ 0.2 pc dirrs would be expected to exhibit a sig- 
nificantly larger binary fraction, which is of the order of 
70 — 80 per cent. 



4-3.4 IGBDFs for galaxy types 

The above results hold if cluster formation in the MW is rep- 
resentative for other galaxy types (Sec. 13.31 rh = 0.2 pc, but 
with different SFRs). Then, the expected period and mass 
ratio IGBDF for ellipticals, spirals and dIrrs for a population 
of binaries with mi ^ 2Mq (no cuts in period or primary- 
mass) are expected to be as in Fig. [O] The period IGBDF 
of Es peaks at shorter periods than the corresponding dis- 
tributions for S and dIrr galaxies. It is due to the initial 
star burst in Es (high SFR), namely that more dynamically 
evolved binary populations from high-mass clusters, which 
are not present in the lower SFR spirals and dIrrs (eq. |6]), 
contribute to the galactic field of Es (Sec. 13. 2p . The mass 
ratio IGBDFs for all morphological types look alike, a dif- 
ference only appearing through the varying binary-fraction 
between the galaxies (the area below the distributions). 



4.4 Model Limitations 

The validity of the results and predictions outlined in this 
analysis depend on the accuracy of the assumptions entering 
the IGBDF model. While the model formulated in Sec. [2:21 
appears robust, the results obtained by performing the in- 
tegration in eq. (|15p depend on the prop erties of the bir th 
and initial binary population (see Sec. [2l lKroupa[[T995a[ lbl. 
Paper I) and the resulting analytical description of the evo- 
lution of binary properties in the A^-body models of Paper 
I, which use this initial binary population. The two ma- 
jor assumptions are that (i) the birth binary population is 
formed via random-pairing of the two binary components 
from the canonical stellar IMF and that (ii) the birth pop- 
ulation evolves into the initial population via pre-main se- 
quence eigenevolutio n, which was sp ecifically parameterised 
for G-dwarf binaries l|Kroupa|[l99"5bh . If, e.g., instead of the 
declining mass ratio dist ribution found in DM9 1 the flat 
mass ratio distribution in [Ragh avan et al.1 (|2010f ) for solar- 
type stars were correct, this would possibly imply a birth 
pairing method for binaries which is different from random- 
pairing and, in turn, wou ld require alteration of the binary 
birth population (see also [Kouwenhoven et al.[[2009[ l. Alter- 
natively, the eigenevolution model might need adjustments, 
e.g. through primary star dependent eigenevolution param- 
eters A = A(mi) and x ~ x('^i)) so that the initial binary 
population may be mildly different for other spectral types. 
These limitations should always be kept in mind, but the 
overall results of the IGBDF model would not be affected. 



5 SUMMARY &: OUTLOOK 

Following observational evidence which implies all stars less 
massive than about 2Mq to form as binaries with compo- 
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nent masses picked randomly from the stellar IMF in dis- 
crete star formation events (i.e. embedded star clusters) 
and allowing for pre-main sequence eigenevolution (Sec. [2|) 
and stimulated evolution, the concept of integrated galactic- 
field binary distribution functions (IGBDFs) is introduced. 
Adding up the stellar populations that ever formed in star 
clusters which are selected from an ECMF yields the statisti- 
cal binary properties of a galactic field population (Dynam- 
ical Population Synthesis). This approach is similar to the 
IGIMF theory which adds up the stellar IMFs in individual 
clusters to calculate the galaxy-wide IMF (Sec.[l]). 

The IGBDFs depend on the minimum cluster mass, 
Afeci.min, the galaxy-widc SFR, the steepness of the ECMF 
and the typical, or average value for cluster half-mass radii, 
rh, in a star cluster system. The galactic field binary- fraction 
increases with decreasing Afcci,min and SFR, with increas- 
ing index f5 of the ECMF (eq. [2|) and with increasing rh 
(Fig. [2]). It is found that low- mass (i.e. low-density) clusters 
contribute most binaries to the field since stimulated evolu- 
tion is least effective in them and low-mass clusters are most 
numerous due to the steep ECMF. High-mass (high-density) 
clusters donate most single stars since binary destruction is 
effective and the number of stars forming in a high-mass 
cluster is larger than in a low-mass cluster. 

Applying the IGBDF model to the MW, i.e. estimat- 
ing Meci.min, P Slid the SFR from observations, the pe- 
riod, mass-ratio and eccentricity IGBDFs are constructed 
for different typical cluster sizes, rh. The models indepen- 
dently agree with observed distributions for late-type bina- 
ries in the solar neighbourhood, solely adjusting the typical 
Th to 0.1 — 0.3 pc, which is the single remaining free model 
parameter. This suggests that MW clusters typically form 
quite compact. M-dwarf binaries appear to have formed in 
slightly more extended clusters than G-dwarfs, a possible 
sign of mass-segregation. The integrated populations show 
that the majority of all Galactic field primaries end up being 
single stars despite being born in binary systems. In partic- 
ular, the Galactic field binary fraction is a function of the 
spectral type, i.e. /b(M) < /i,(K) < /b(G) < fb(F), and the 
binary-frequencies derived from binary-star formation and 
stimulated evolution in compact clusters agree with obser- 
vational data. 

It has been pointed out that the shape of the mass 
ratio distribution depends on the considered primary-mass 
(Fig. [5]). While F- to K-type binaries show a decreasing trend 
with increasing q, according to the model the M-dwarf mass 
ratio distribution increases as q increases. A fiat distribu- 
tion is only obtained if primaries from the full mass-range 
are combined when constructing the mass ratio distribution. 

Using the best-fitting model, the named integrated dis- 
tributions for late-type binary populations (F to M) in 
the solar neighbourhood are predicted (Fig. |4|. While the 
IGBDFs for F-, G- and K-type binaries are probably hard 
to distinguish by observation, M-dwarf binary- and the all- 
binary cumulative distributions appear rather distinct from 
them. Therefore the M-dwarf and cumulative distributions 
are probably the best populations, in comparison with the 
ones for F-, G- and K-binaries, to test the IGBDF predic- 
tions. 

Assuming star formation in the MW (in terms of the 
ECMF and cluster radii) is typical also for other galaxy 
morphologies, the binary population in elliptical galaxies 



(high SFR) is predicted to be significantly smaller than in 
spiral (intermediate SFR) and dwarf galaxies (low SFR) and 
their period- and mass-ratio IGBDFs are calculated. 

The IGBDF model will be extended to allow calcula- 
tion of binary populations in galaxies which had a strongly 
varying star formation history, i.e. a time-dependence 
will be incorporated, via SFR(t). It will be possible to 
synthesize binary properties in individual nearby (dwarf) 
galaxies, and the case of the MW may be revisited. A full 
synthetic galaxy might be constructed and "observed" in 
the computer to mimic real data in order to provide more 
sophisticated means of comparison with observations. It 
will also allow to test their influence on observationally 
derived properties, especially in populations where binaries 
cannot be resolved, such as velocity dispersions (used to 
calculate dynamical masses). 
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